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10-Nov-85 



Subject: Mandelbrot Sets 



From: Russ Altman 



I have uploaded the basic program for native /// mode that will produce 
Mandelbrot Sets in color 190 * 140 mode. The description reminds you to do 
>perform grafixon before running the program, so you can watch its progress. 
It takes on the order of an hour + to finish a picture, so it is a very good 
idea to turn off the display refresh (CONT-KEYPAD-5) to speed up the execution 
(you can put it back on every now and then to admire the picture as it is 
built) . 

If it does not work, please let me know SOON, since it worked before I 
uploaded it, and I have used it a number of times. However, I must admit 
that I added some comments from a text editor, and if I did something stupid, 
there might be problems. 

A key thing to realize about the program is that by changing the parameters 
xhi, xlo, yhi, ylo, you can "zoom in" on sections of the image that may be 
interesting. There are a lot of interesting patterns at the border which 
deserve close scrutiny. Enough said, let me know how it goes... Russ 
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* MANDELBROT SET PLOTTING PROGRAM FOR APPLE /// ******* 
Program requires "BGRAF . INV" takes 1-2 hours to complete 
an image (See Scientific American, October 1985) . 
Please forgive bad basic form, I don't usually program 
in this language. Russ B. Altman 11/5/85 

XHI, YLO, YHI are boundaries of the box to be drawn on screen 

starting values are -2,2 for both sets. Can change to 
closeups of certain areas. 

.Dl/BGRAF . INV" :REM Change this path depending on where you keep 
your Basic invokables . 



REM **** 
REM 
REM 
REM 
REM 

REM XLO, 
REM Best 
REM get 
INVOKE 
REM 
XLO=-2 
XHI=2 
YLO=-2 
YHI=2 

REM XMAX = Maximum number of points in X direction 
REM YMAX = Maximum number of points in Y direction 
XMAX=139 
YMAX=191 

REM XMAX and YMAX above are resolution of screen 
REM Now compute the sample size DELTAX, DELTAY 
DELTAX= (XHI-XLO) /XMAX 
DELTAY= (YHI -YLO) /YMAX 
PRINT DELTAX, DELTAY 

REM Set max to be the maximum number of colors you can draw at each dot 

MAX%=15 

hits%=0 

REM Now call all the functions to initialize graphics, clear screen 
REM the ones below are for Apple /// 

PERFORM INITGRAFIX 

PERFORM GRAF IXMODE ( % 3 , % 1 ) 

PERFORM fillcolor(%0) 

PERFORM FILLPORT 

REM Now just loop through all values of I, J and calculate corresponding 
REM values of X, and Y, then GOTO the subroutine that calculates ANS . 
REM Depending on what ANS is, either draw a dot of color ANS or not 
FOR I%=1 TO INT (XMAX) 

PRINT"I = ",I% 

FOR J%=1 TO INT (YMAX) 
GOSUB 300 

REM If ANS not equal to 0, then draw a dot at the point I, J 
IF ANSO0 THEN GOSUB 230 
ANS=0 
NEXT J% 
NEXT 1% 
STOP 

REM This is the first subroutine which just draws a dot at 
REM point I, J of color ANS . make sure arguments to functions 
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REM are right data type (string, int, real, etc..) 

XINT%=I% 

YINT%=J% 

REM Make sure we have an integer for color 
A%=INT(ANS) 

REM PENCOLOR sets the color of the dot to be drawn (on Apple ///) 

REM DOTAT actually draws a dot at XINT, YINT = I, J 

PERFORM PENCOLOR (%A%) 

PERFORM DOTAT (%XINT%, %YINT%) 

HITS%=HITS%+1 

REM HITS just keeps track of how many non-zero dots — expendable 
RETURN 

REM this subroutine actually calculates ANS, or how long it takes 

REM for a number to be squared and reach a limit (4) 

X=I%*DELTAX+XLO 

Y=J%*DELTAY+YLO 

FIRSTX=X 

FIRSTY=Y 

REM N is a counter for how many times we have calc'ed MOD2 

N=0 

N=N+1 

IF N>MAX% GOTO 4 40 

REM If MOD2 > 4 then we are done, so leave subroutine 

REM else, we want to calculate a new set of X,Y'S and recalc MOD2 

MOD2=X*X+Y*Y 

IF MOD2>=4 GOTO 4 60 

oldx=x 

oldy=y 

x=oldx*oldx-oldy*oldy+f irstx 
y=2*oldx*oldy+f irsty 

REM Loop back up and try calcing MOD2 again to see if reached limit 
GOTO 350 

REM If we haven't reached MOD2 by the time we have used up all 
REM possible colors, then we just draw black (return ANS = 0) 
ANS=0 
RETURN 

REM We set ANS to N, so that we know what color to draw. Note that 

REM N and ANS refer to the number of times we calced MOD2 *and* 

REM the color we are going to use to plot on the screen. 

ANS=N 

RETURN 

REM End of program Mandelbrot .Bas 
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